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Abstract 

Using as an illustrative example the p=l operator-ordered Wheeler-DeWitt equation for 
a closed, radiation-filled Fricdmann-Robertson- Walker universe, we introduce and discuss the 
supersymmetric double Darboux method in quantum cosmology. A one-parameter family 
of "quantum" universes and the corresponding "wavefunctions of the universe" for this case 
are presented. 
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Many interesting results have been obtained mostly in one-dimensional quantum mechanics 
by means of Darboux- Witten (DW) supersymmetric procedures; for a recent review see 
These are, essentially, factorizations of the one-dimensional Schrodinger differential operator, 
first performed in the supersymmetric context by Witten in 1981 [||, and in mathematics lit- 
erature, in the broader sense of the Darboux transformation, by Darboux, as early as 1882 

Constructing families of isospectral potentials is an important supersymmetric topic which may 
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have many physical applications Q. Since in quantum cosmology the Wheeler-DeWitt (WDW) 
equation is a Schrodinger-type equation one might think of introducing the isospectral issue 
in this field as well. This is our purpose in the present work. To accomplish it, we shall use 
the simple example of a radiation-filled, closed Friedmann- Robertson- Walker (FRW) "quantum" 



universe. 



We briefly review the well-established supersymmetric double DW technique of deleting fol- 
lowed by "reinstating" a nodeless, radial (half-line) ground-state, say uo(a), of a potential V~(a) 
by means of which one can generate a one-parameter family of isospectral potentials Vi so {a\ A), 
where A is a labeling, real parameter of each member potential in the set. The family reads 
explicitly 



V lso (a; A) = V-(a) - 2[ln(J + A)]" = V~(a) - + ^ , (1) 

where the primes are derivatives with respect to the radial variable a and 

Jo(a)= [ a u 2 (y)dy. (2) 
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The result of Eq. (1) is obtained if one factorizes the one-dimensional Schrodinger equation 
with the operators A = -j| + W(d) and = — + W(a), where the superpotential function 
is given by W = —Uq/uq. Then, in the DW scheme the potential and superpotential enter an 
initial 'bosonic' Riccati equation V~ = W 2 — W'. At the same time, one can build a 'fermionic' 
Riccati equation V + = W 2 +W , corresponding to a 'fermionic' Schrodinger equation for which 
the operators A and A + are applied in reversed order. Thus, the 'fermionic' potential is found 
to be V + = V~{a) — 2 (t^J ■ This potential does not have uq as the ground state eigenfunction. 
However, it is possible to "reintroduce" the uo solution into the spectrum by means of the general 
superpotential solution of the fermionic Riccati equation. The general Riccati solution reads 

W gen = W(a) + ^ ln[J (a) + A] . (3) 

The way to obtain Eq. (3) is simple and well-known ||. From Eq. (3) one can easily get Eq. (1) 
by just using the bosonic Riccati equation. But, is there any price for this procedure ? Yes, there 
is. The new wavefunction differs from the initial one by a 'normalization', A-dependent constant. 



The important thing which was noticed in the literature was the damping nature of the integral 
Jo, both for the family of potentials as for the wavefunctions. The parameter A is acting as 
a "damping distance" which is signaling the importance of the integral term. Finally another 
useful formula gives the family of A-dependent (non-normalizable) wavefunctions as follows 

u (5) 



(a; A) 



(4) 



Jo + A 

We choose now one of the most simple cosmological metrics, a closed FRW model filled only 
with radiation, in order to clearly illustrate the supersymmetric isospectral issue in quantum 
cosmology. Kung Q has recently discussed the WDW equation for the closed, k > FRW 
metric with various combinations of cosmological constant and matter. The WDW equation 
used by Kung reads 
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u(a) = , 



(5) 



where the variable a is the cosmological scale factor. The parameter p enters as a consequence of 
the ambiguity in the ordering of a and d/da; p = 1 is the so-called Laplacian ordering, whereas 
p = 2 is convenient for the WKB approximation. 

In the matter sector various contributions may be present but we shall consider only the 
radiation case, i.e., 
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where p r is the radiation energy density. The most compact form of the WDW equation in the 
FRW metric can be written down as follows, 



d _ d 



+ V(a) 



u(a) = , 



(7) 



da da 

where the tilde variables are rescaled and dimensionless ones. Thus, only with the radiation 
matter sector the cosmological potential turns out to be an oscillator potential of the form 



V(a) = a 2 - P 2 , 
3 



(8) 



where a 2 = ^j^a 2 is the tilde scale factor of the universe, and 2 = ^-p r dQ expresses the 
radiation effect on the cosmological expansion. 

With the ansatz u(a) = g(a)e~ a2 1 2 , and with the further change of variable x = a 2 , one gets 
the following confluent hypergeometric equation for g(x), 



Using the p = 1 factor ordering and with the natural assumption that the physical u must 
vanish asymptotically for x —* do, Eq. (9) turns into an eigenvalue problem for the Laguerre 
polynomials g{x) = L° n _^- 2 _ 2 y^(x), where n G N + . Thus, in this case, the wavefunctions of the 
universe are of the type 



u n {a) = e-' a2 ' 2 L^ 2 _ 2)/A {a 2 ) . (10) 

One can also see that the radiation energy density is quantized through the condition p 2 = 4n+2. 
Thus, = y/2 is the lowest value, corresponding to the Gaussian ground state wavefunction 
wo (a) = e _fl2//2 . Other factor orderings lead to various confluent hypergeometric functions g and 
to a change of the quantization condition [|| . 

In the following, we shall perform the double Darboux construction using the Gaussian node- 
less wavefunction of the closed, radiation-filled FRW universe. One can easily obtain W = a, 
that is the superpotential is the dimensionless scale factor of the universe, and V + = V~ + 2 = 
a2-p2 + 2 = a 2 . 

In this case, the one-parameter family of bosonic cosmological potentials can be written as 
follows, 

V iso = a 2 -2+ f ,-x° + n r x , r 11 
A s erf (a) + A (A s erf (a) + X) z 

In order to avoid singularities one should take |A| > X s = |^|. We have plotted some 

members of the family equation (11) in Fig. 1. For A — ► ±oo the members are very close in 

shape to the original V~ potential. 
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The one-parameter family of wavefunctions of the universe are 

Ugen(a; A) = - — y° . (12) 
A s eri(a) + A 

The members of the family of wavefunctions given by Eq. (12) are plotted in Fig. 2 for the same 
A values as for the potentials. They can also be written as "Gaussians" u gen = e -°- 2 l 2rj2 ( a ) with 
a dispersion depending on the scale factor of the universe 

a 2 

a2{d) = 5 2 + ln(A s erf(a) + A) 2 ' (13) 
We recall that Kiefer || has considered wavepackets with time-dependent dispersion within the 
WKB approximation to minisuperspace models. 

In conclusion, we have introduced concepts like isospectral cosmological potentials and a 
corresponding one-parameter family of wavefunctions of the universe for the simple case of 
closed, radiation-filled FRW "quantum" universe. 

The work must be considered as merely a brief introduction of the isospectral supersymmetric 
problem in quantum cosmology. Its further applications in this field will be fully explored in 
future investigations. 

This work was partially supported by the CONACyT Project 4868-E9406. 
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Fig. 1 

Members of the one-parameter isospectral family of closed, radiation-filled FRW potentials for the 
following values of the A-parameter (from the right to the left): (a) -.8863; (b) -.8903; (c) -.9303 
together with the original potential (the evolvent) V~— = a 2 — 2. 
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Fig. 2 

Squares of the "wavefunction of the universe" for the same values (right to left) of the A parameter as in 

Fig. 1. The largest peak is scaled down forty times. 
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